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We consider controlled collisions between two ultracold atoms guided by external harmonic po- 
tentials. We derive analytical solutions of the Schrodinger equation for this system, and investigate 
the properties of eigenergies and eigenstates for different trap geometries as a function of a trap 
separation and of the scattering length. When varying the trap separation the energy spectrum 
exhibits avoided crossings, corresponding to trap-induced shape resonances. Introducing an energy- 
dependent scattering length we investigate the behavior of the system in the vicinity of a magnetic 
Feshbach resonance. Finally, we illustrate our analytical results with two examples: the quantum 
phase gate controlled by the external magnetic field, and a scheme for a coherent transport of atoms 
in optical lattices into higher Bloch bands. 
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I. INTRODUCTION 



Atomic interactions in the presence of an external con- 
finement represent one of the major ingredients for the 
schemes implementing quantum information processing 
in systems of neutral trapped atoms (see e.g. m for a re- 
view). The tight external confinement for neutral atoms 
can be created by means of optical lattices Q , atom chips 
[1, Q and dipole traps Q . Realization of Mott insulator 
phase 0, 0, S| gives the possibility to precisely control a 
number of atoms confined in a single well of an optical 
lattice. Moreover, magnetic Feshbach resonances, widely 
used in recent experiments on ultracold atoms [9[, per- 
mit for arbitrary tuning of atomic interactions, which was 
the key ingredient to accomplish molecular Bosc-Einstcin 
condensates [loj and the superfluidity in ultracold Fermi 
gases All these achievements make systems of ul- 

tracold neutral atoms very attractive in the context of 
quantum information processing and quantum state con- 
trol. This has been demonstrated in experiments on cre- 
ation of massive entanglement [l^ . realization of quan- 
tum gates or a coherent control of atoms in op- 
tical lattices. Recent studies of cold ion-atom collisions 
[isl . [TgI opens also a way for the future realization of 
quantum computation schemes involving both atoms and 
ions 113 J combining the advantages of strong inter- 
actions of charged ions with long decoherence times of 
neutral particles. 

Modeling of ultracold atomic interactions is most con- 
veniently done in terms of an s-wave delta pseudopoten- 
tial [l^, [l^ . In this approach the system of two inter- 
acting atoms in a harmonic trap has analytical solutions 
for spherically symmetric pl| or axially symmetric har- 
monic potentials [l^, HI]- Moreover, generalization of 
the pseudopotential to higher partial waves [l^, [H, [l^ 
allows to solve the problem for generic types of central 
interactions in the presence of the harmonic confinement 
[H, m, [13, [13 ■ The system properties become even 



more intriguing when the two particles are trapped in 
separate harmonic potentials. In such a case the sys- 
tem exhibit trap-induced shape resonances [2^ between 
molecular and trap states. These resonances, manifesting 
itself as avoided crossings in an energy spectrum, can be 
applied, for instance, for a quantum state control [13, [13; 
or realization of quantum gates [13 . 

The setup consisting of two atoms in separate traps can 
be realized with spin-dependent lattice potentials [soj . 
allowing for individual control of the trapping potential 
for neutral atoms depending on their internal hyperfine 
state. Another example is the system of a single atom 
and a single ion that can be trapped in independent 
potentials created by optical-dipole and electric radio- 
frequency fields [l3|. In the latter case the long-range 
nature of the ion-atom potential prevents from using 
the contact pscudopotentials to model the interactions, 
nevertheless, supplementing the pseudopotential with an 
energy-dependent scattering length [3l|, [13, allows to 
consider the regime when the external confinement is of 
the order of the interaction range. For neutral atoms 
application of an energy dependent scattering length ex- 
tends the validity of the pseudopotential treatment to 
the case of very tight traps or large scattering lengths in 
the vicinity of resonances [H, [13 . It also properly ac- 
counts for the whole molecular spectrum [20| . Discussed 
models based on the use of pseudopotential provides for 
a very accurate description of neutral atom interactions 
in the presence of external confinement. This has been 
confirmed, for instance, in the recent experiments on the 
creation of homonuclear [3^ and heteronuclear 35 1 dimer 
molecules in optical lattices. 

In this paper we present the exact analytical solu- 
tions for two interacting atoms confined in separate har- 
monic potentials. So far the eigenenergies and eigen- 
states of such a system have been studied only numeri- 
cally [13, [13, [13 ■ For simplicity we assume that atoms 
are confined in the traps of the same trapping frequency, 
which allows to separate the center-of-mass and rela- 
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FIG. 1: (Color online) a) Two atoms in cylindrical harmonic 
traps separated by the distance d. b) A relative part of the 
effective potential. 



tive motions. The numerical studies of the more general 
situation with different trapping frequencies [l7| show, 
however, that the basic features consisting in the pres- 
ence of the trap-induced resonances remain unchanged. 
We present the analytical results for the energy spec- 
trum and wave functions discussing different geometries 
of the harmonic trapping potential. Applying the energy- 
dependent scattering length we investigate the effects of 
Feshbach resonances on the trap-induced resonances. We 
illustrate our analysis with two example applications of 
the trap-induced shape resonances for the atoms in op- 
tical lattices. We consider a simple scheme for a quan- 
tum phase gate between atoms in separate wells, and we 
present a method for a coherent transport of atoms into 
higher Bloch bands. 

The paper is organized as follows. In section |TT] we de- 
rive analytical solutions of the Schrodinger equation for 
two interacting atoms in separate traps. Scction lTTTl is de- 
voted to the analysis of the energy spectrum. The wave 
functions arc investigated in section llVl The energy spec- 
trum in the vicinity of a magnetic Feshbach resonance 
is analyzed in section Ifff Dl We illustrate the applica- 
bility of our analytical results in section |V] considering 
the quantum phase gate and the scheme for a coherent 
transfer between quantum states in the trap. Section IVll 
present some conclusions, and three appendices give some 
more technical details on our derivation. 



II. MODEL 

We consider the system of two interacting particles 
confined in separate harmonic potentials. The setup is 
illustrated schematically in Fig. [TJ Wc assume that the 
trapping potentials of the particles can be controlled in- 
dependently and the distance between traps can be arbi- 
trarily adjusted. The Hamiltonian of the system can be 
written down as: 



H 



2mi 2to2 
+ F/''(ri - dl) + V^''\y2 - d2) H- F,(ri - V2). 



(1) 

where mi, are masses of atoms, v/"'^''(r), V^'^'(r) are 
their respective trapping potentials, di, d2 denote trap 



positions, and T^(r) is the interaction between particles. 
For simplicity we assume that the harmonic traps are 
cylindrically symmetric, however, our method to derive 
analytical solutions is general and can be applied for ar- 
bitrary harmonic potential. The trapping potential can 
be written as 



.7-1,2. 



(2) 



In a harmonic external potential the centcr-of-mass and 
relative motions can be separated, and in the following 
we will focus only on the relative motion that is governed 
by 



Hrel = 



1 



2/X 



"sj^ + -lJ.uj\v~dY + V{v). 



(3) 



Here, r = ri — r2, denotes the relative coordinate, ji = 
mim2/{mi +m2) is the reduced mass, and d = di — d2 is 
the distance between the traps. Without loosing gener- 
ality we can assume that traps are displaced along z- 
direction d = dz. In our model we approximate the 
atom-atom interaction with a contact Fermi pseudopo- 
tential 



l/(r) 



/i or 



(4) 



describing s-wave scattering, which dominates at ultra- 
cold temperatures. This approximation is valid provided 
that the characteristic length of the trapping potential 
is much larger than the characteristic radius of the in- 
teraction. This is typically well fulfilled for neutral atom 
collisions, however it typically breaks down for long-range 
potentials, like the atom-ion one. In such cases one can 
apply the energy-dependent scattering length to extend 
the validity of the pseudopotential to the regime when 
the trap size is of the order of the interaction range. 

In the following wc will adopt dimensionless units ex- 
pressed in terms of the oscillator length I ~ ^Jhllpmoz) 
and energy hio. Rewriting the Hamiltonian in cylindrical 
coordinates yields 



-iA + i(/7V' + df ) + 27ra(S(r)^r 



where rj = to x^juj^ is the trap aspect ratio. 



^'(r) = £;*(r), 
(5) 



III. ENERGY SPECTRUM 

We turn now to evaluation of eigencnergics. To this 
end we decompose Vl'(r) in the basis of the harmonic os- 
cillator wave functions in cylindrical coordinates. Since 
the delta pseudopotential acts only on states that do not 
vanish at r = 0, in the expansion we include only the 
harmonic oscillator states with zero projection of the an- 
gular momentum on the symmetry axis (m = 0). The 
expansion reads 



(6) 
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where $„.o(PiV) ^-^d &kiz) are the harmonic oscilla- 
tor states for one and two dimensions, respectively, in 
the trapping potential Vt(r). Substituting ([6|) into the 
Schrodinger equation ([5]), and projecting on a single state 
of the basis gives the expansion coefficients 



c 



^;,o(o,^)efeM) 

E — Enk ' 



(7) 



where En^k = 1/2 + 77 + ^-1- 27?n are eigenenergies for a 
cylindrically symmetric trap and 



C = 27ra 



d 



— (r ^ Cn',k,^7i'fl{p, ip),Qkiz - d)) 



(8) 

is some constant prefactor, given by the normalization 
of the wave function. Substituting Eq. ([7]) into ([8|) wc 
obtain equation determining the energy spectrum 



1 

27ra 



d 



rvl/g(r) 



(9) 



where \f e (r) is an eigenstatc corresponding to the eigen- 
energy e 



$:,o(0, <^)$n.o(p, y)9^(-d)9fc(z - d) 
27]n + k — e 



(10) 

Here, e — E — Eq and Eq = 1/2 + i] is the energy of the 
zero-point oscillations. Wc note that the eigenstatc (fTO|) 
is not normalized. We further proceed by expressing the 
denominator in Eq. (jlOp in terms of an integral 



1 



27771 + k + e 



(11) 



which allows to perform summations over excitation 
quantum numbers k and n. To this end we substitute 
the explicit formulas for the oscillator wave functions 
$„_o(/3, (f) and 6fc(z), and utilise the summation formulas 
for the Hermite and Laguerre polynomials (see Appendix 
\X\ for details) . This yields an integral representation of 
the wave function 



(2^)! 



dt- 



■ coth(t)- 



■ coth{r}t) + 



d(d-z) ■, 



■\/sinht sinh(77i) 



(12) 



Similar integral representation can be derived for the 
equation determining the energy spectrum. To this end 
we substitute ([12]) into ([9]), obtaining 



a V 2/ 



dr 



(13) 



where 



-t/2 



T{x) 



dt{ 



— rf " 

rje e 



(l-e-*'')yr 



1 , 

^3/2' 



(14) 



The term t~^/^ assures that the integral converges at the 
lower limit of integration, and it results from the appli- 
cation of the regularization operator in the Fermi pseu- 
dopotential. As can be easily shown, the wave function 
(|10p exhibits 1/r divergent behavior at small r, origi- 
nating from the behavior of the integral at t — > 0. 
Factoring out the divergent term from the integral rep- 
resentation p2)) in a way described in Ref. [2^ , and then 
removing the divergent 1/r term with the help of the reg- 
ularization operator results in the integral In the 
simplest case of an isotropic trapping potential (77 — 1), 



the integral (fH)) can be evaluated in terms of a series 



J^{x,r] = 1) 



-2V^. 



r(x) 



r(x-i) 



+J2 ^^i!^^(2.T, A;-i)2Fi(2x-2, fc-i; 2x+k^^; -1), 

k=l 

(15) 

where B{x,y) = T{x)r{y)/T{x + y) denotes the Euler 
Beta function [H]. For d = this straightforwardly re- 
duces to the known result for two atoms in isotropic trap 
[2H . In general case the integral (|14p is too complicated 
for an analytical treatment. While for x > the inte- 
gral is well defined and can be calculated numerically, 
for a; < the integral is defined only in the sense of 
the analytic continuation. Therefore, in the following we 
will develop the recursion formula, which relates T{x) at 
2; < with the values of J^{x) at x > 0. 



A. Recurrence relations 

Let's consider the difference between values of lP{x) 
separated by ry 

^2 1-e-*/2 



T[x) - T(x + r/) 



dt 



(16) 
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The above integral can be calculated analytically (see 
Appendix B for details) 

Tix)~-T{x + T]) = B (i, 2a;) $i (i, 2x, 2x + i; -1, -d^) , 

(17) 

where $i(a, 6, c; z, w) is the degenerate hypergeometric 
function of the two variables. It is defined as 13611 



^i{a,b,c;z,w) = ^ 



71 



(18) 



where (a)„ is the Pochhammer symbol. Although the in- 
tegral representation (fT4|) is valid for a; > 0, we note that 
the final result involving $i(a, 6, c; z, w) is well defined 
for both a; < and a; > 0, by the virtue of the analytic 
continuation. While $i allows to express the result of in- 
tegration in quite elegant form, its numerical evaluation 
is rather cumbersome. In fact, in the numerical calcula- 
tions we have applied the following series expansions 

B (i, 2x) $1 (i, 2x, 2a; + i; -1, -d^) 

J{n + 2x) 



E 

n=0 



(-ir 



r(2x)n! 



X B(2a;,7i+ i)iFi(n+ i;2a; 



r(2a;) ^ n!r(2n + 4a;) 4' 



which are derived in the Appendix [ 



i- -d^ 

9 7 " 



1 

2' 



(19) 
(20) 



B. Results 

We calculate the energy spectrum from Eq. ([13]) with 
T(x) evaluated numerically from the integral (|14p for 
a; > and using the recurrence relation (jl7p for a; < 0. 
Fig. [5] shows the results for the spherically symmetric 
traps and a fixed value of the scattering length, but with 
varying the trap separation d. We observe that for large 
d the energy levels acquire values of a noninteracting har- 
monic oscillator. For smaller trap separations we observe 
avoided crossings with a bound state of the interparti- 
cle potential, which is lifted up by the external trapping 
potential, with approximate quadratic dependence on d. 
This behavior can be explained by observing that the 
bound state wave function is concentrated around r = 0, 

so {^raoliMHreMWraolWl ^ Eb + ^flLU^^d^, whcrC Eb 

denotes the binding energy. 

The avoided crossings correspond to the trap-induced 
shape resonances [2^, that occur when the energy of 
bound states lifted by the trapping potential coincide 
with the energy of a vibrational state in the trap (see 
Figfflb). For the spherically symmetric traps (77 = 1) at 
zero separation {d = 0) the angular momentum of the rel- 
ative motion is conserved and the states can be addition- 
ally characterized by the partial wave quantum number I. 
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FIG. 2: (Color online) Energy spectrum (solid red lines) 
versus trap separation d for two interacting atoms confined 
in harmonic traps of geometry 77 — uj±/uj2 = 1, and for 
a w 0.79/. All lengths are expressed in oscillator units; 
I — ^JhJijjuJz)- Blue dashed line shows the approximate 
quadratic shift of the bound state due to the external trapping 
potential. Black dotted line is a result of a perturbation cal- 
culations between exact spherical zero-separation wave func- 
tions and oscillator p states. 



In Fig. [5] we observe that at d = the spectrum contains 
not only the s-wave eigenstates with energies shifted by 
the zero-range potential plj . but also p-wave harmonic 
oscillator states with energies En — fiio{h/2 4- 2n). The 
presence of the p-wave states in the spectrum can be ex- 
plained in the framework of the perturbation theory. At 
small d the external potentials i/xw^(a;^ + + (z — d)^) 
couples s-wave eigenstates of the interacting atoms with 
p-wave harmonic oscillator states (see Fig. [2]). 

Similar features at small d can be also observed for 
an anisotropic trap (see Fig. [3]). In this case eigen- 
states of the s-wave contact potential are coupled in the 
first order in d to the harmonic oscillator states. In the 
cylindrically symmetric trap their energies are given by 
E = hhj{ri + 1/2 + Uz + 2r]n±). By studying the case of 
non-commensurable trapping frequencies tOz and lo±, we 
have verified that only states with even Uz are present in 
the spectrum at d ^ 0. 

Figs. [3] and [5] present the energy spectrum versus the 
scattering length, for a fixed trap separation, and for two 
different trap geometries: rj = 2 and rj ~ 0.7, respec- 
tively. For d — the spectrum coincides with the one ob- 
tained for two trapped atoms in an anisotropic harmonic 
trap [23j . For non-commensurable trapping frequencies 
T] = 0.7 and d = 0, we note the appearance of some ad- 
ditional states in the spectrum, weakly dependent on a. 
These states give rise to avoided crossings at a ^ (23t . 
For small and positive d the energy spectrum contains ad- 
ditionally states of the harmonic oscillator with even , 
which may be observed as horizontal lines, due to their 
weak dependence on a. They create avoided crossings 
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FIG. 3: (Color online) Energy spectrum (solid red lines) 
versus trap separation d for two interacting atoms confined 
in harmonic traps of geometry 77 — uj±/uJz = 2, and for 
a ~ 0.79?. All lengths are expressed in oscillator units; 
I — ^/h/{ij.uiz). Blue dashed line shows the approximate 
quadratic shift of the bound state due to the external trapping 
potential 
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FIG. 5: (Color online) Energy spectrum as a function of the 
inverse scattering length for two interacting atoms confined 
in harmonic traps of geometry rj — uj±/uJz = 0.7, and for trap 
separations d/l — 0.0,0.1,2.1,4.0. All lengths are expressed 
in oscillator units: / = ^/KJJpuj^. 



C. Energy-dependent scattering length 
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FIG. 4: (Color online) Energy spectrum as a function of the 
inverse scattering length for two interacting atoms confined 
in harmonic traps of geometry rj — ui±/io^ — 2, and for trap 
separations d/l = 0.0,0.1,2.1,4.0. All lengths are expressed 
in oscillator units: I — \/h/ (hlJz). 



with the rest of the states that are strongly influenced by 
the presence of the atom-atom interaction, and which are 
already present at c? = 0. With increasing d, the energy 
of the bound state rises, and at some moderate trap sep- 
aration (d ~ the bound state starts to create avoided 
crossings with vibrational states in the trap. 



So far we have performed our analysis using the stan- 
dard, energy-independent scattering length, defined in 
the zero-energy limit: a = — lim;;_to tan 6Q{k) / k in terms 
of the s-wave phase shift So{k) and the wave vector 
k = \/2ilE/h. This approach is valid for sufficiently 
weak trapping potentials: fuu <C / {ji\a\Rcs) . In this 
case the effective range correction in the expansion of 
the s-wave phase shift can be neglected: kcoi5[j{k) = 
-l/a + ii?offfc^ + 0(fc^), where we can take ?i^fc^/(2/x) ~ 
hw. Here i?off denote the effective range, which for a 
van der Waals interaction is of the order of the charac- 
teristic length of the C^/r^ potential: _Rcff ^ Rvdw = 

i (2^^6/^2)1/4 jia. 

In the situation when the ccnter-of-mass and rela- 
tive motions are not coupled, the more general treat- 
ment of atom-atom interactions is based on the use of an 
energy-dependent scattering length a{k) = —taiiSo{k)/k 
jsil . [32 | . This is particularly important for large values of 
the scattering length obtained in the presence of scatter- 
ing resonances, that may occur due to the interchannel 
coupling (Feshbach resonances) or due to the potential 
barrier (shape resonances). Moreover, the energy depen- 
dent scattering length extends the validity of the pseu- 
dopotential to the case of tight traps. In such system the 
standard, energy-independent pseudopotential becomes 
inaccurate because of the relatively high kinetic energy 
of the colliding atoms and comparable range of the inter- 
action potential and of the external confinement. 

We illustrate this behavior in Fig. [HI comparing pre- 
dictions of the energy-independent and energy-dependent 
pseudopotentials for a model square-well potential in- 
teraction. We note that for a narrow well of the size 
Rq = 0.05/, the standard contact potential is quite ac- 
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FIG. 6: (Color online) Upper panels: energy spectrum ver- 
sus trap separation d for two interacting atoms confined in 
spherically symmetric trap for a 0.79L Bottom panels: en- 
ergy spectrum versus inverse scattering length 1/a for atoms 
in spherically symmetric traps separated by distance d = l.ll . 
All lengths are expressed in oscillator units: I = yi'h/{ij,ui^). 
Black solid curves show the result for the energy-independent 
pseudopotentials. Green (light grey) solid lines show the en- 
ergy levels values for the energy-dependent pseudopotential 
with the phase shift given by the square well potential. The 
left panels are calculated for a square well with diameter 
Rf) = 0.05/ and depth V = 520 hui, while the right panels 
are obtained for _Ro = 0.5/ and V = 520ftw. 



curate, with exception of a deep bound state shown in 
the panel c). In contrast, for a wide square well with 
Rq = 0.5/ the two approaches totally differ with respect 
to the energy dependence of the bound state, and as a 
consequence the positions of the trap-induced shape reso- 
nances shown in panel b). The differences are also visible 
in the panel d) in the vicinity of the resonance (a = ±oo), 
where the standard pseudopotential breaks down due to 
the lack of the effective range correction. 



D. Feshbach resonances 

Applying the concept of the energy-dependent scatter- 
ing length, one can analyze behavior of the system in the 
presence of Feshbach resonances. To this end we utilize 
the following result for the energy-dependent scattering 
leng th in the vicinity of a magnetic Feshbach resonance 

m 



[id, l38| , starting from the Breit-Wigner formula for the 
phase shift near a resonance Eq. (PT|) accounts for 
the effects of the finite kinetic energy of colliding atoms 
and applies both to the entrance-channel dominated (so 
called "wide") and closed-channel dominated (so called 
"narrow" ) Feshbach resonances. The main correction due 
to the finite kinetic energy in the open channel is a shift 
of the resonance SB = E/fim, which is included in the 
denominator of Eq. ((2T|) . The term ABE/Eb represent 
the correction, which, as we show below, is important 
for the case of a narrow resonance, whereas it can be 
safely neglected in the case of a wide resonance. Finally 
(1 -I- E/ Eb) in the numerator of ([2T|) describes the mod- 
ification of the resonance width at finite energies. 

A useful parameter that can be used to classify the 
resonances is 77 = {a/ahg)fi'^/{2fia^ABiirn) where a 
is the mean scattering length a = 0.956i?vdw of the van 
der Waals potential. For entrance- and closed-channel 
dominated resonances 77 1 and 77 ^ 1, respectively. 
Since Obg ~ a, the above conditions can be rewritten as 
AB^ni 3> Eb and ABfi,n <C Ei,, respectively ^. Hence, 
for the open-channel dominated resonance one can ne- 
glect the term ABE/Eb with respect to E/^m in the 
denominator of ([21]) . while for the closed-channel domi- 
nated resonances the contribution from ABE/ Eb is im- 
portant. 

For weak trapping potentials: hw <^ Eb and huj <^ 
fimAB, one can omit the corrections due to the finite 
kinetic energy, and in this case Eq. pT|) reduces to the 
well-known result acs{B) — abg [1 — AB/{B — Bq)] for 
the scattering length in the vicinity of a magnetic Fesh- 
bach resonance. In general case, however, one can apply 
the energy-dependent scattering length aefi(£^, -B), which 
properly accounts for the effects of tight trapping poten- 
tials [la, [Si]. 

Here we present some example calculation for a spe- 
cific system of two ^"K atoms. For s-wave collisions 
between hyperfine states \F — 9/2, Alp ~ —9/2) and 
F = 9/2, Alp = —5/2) we use the following parameters 
i^: abg = 174aBohr, AB = 0.97 mT, Bq = 22.421 mT, 
fi^/h = 21.8 MHz/mT, and uj = 2n x 100 kHz. For this 
choice, 77 = 0.4 and this resonance can be classified as an 
intermediate between the entrance- and closed-channel 
dominated limits. Hence, in our calculations we apply 
the full expression ((2T|) for the energy-dependent scatter- 
ing length. 

Figures [7] and [8] show the energy spectrum for the trap 
aspect ratios 77 = 2 and 77 = 0.7, respectively, and for 



aes{E, B) = Obg 



AS(1 + E/Eb 



B-{Bo + E/^im - ABE/Eb 



(21) 

Here, Obg is the background scattering length, Eb = 
/i^/(2/ia^g), is a difference of magnetic moments be- 
tween open and closed channels, Bq is a resonance posi- 
tion and AB denotes its width. This expression can be 
derived from the analytic theory of Feshbach resonances 



The fraction Z of the molecular wave function in the closed 
channel can be obtained from the binding energy of the dressed 
molecular state Z = [1/ fim)dE^al{B) /dB llOll . by applying the 
formula i?niol(^) = / {'^^'■''■eS (^mob ^)^) fo^ the energy of the 
molecule. This yields Z = 1/ [l + |(aefr(S)/abg)A_B/im/_Ei,] , 
which results in the conditions ABfim 2> and ABfj,m "C 
for entrance- and closed-channel dominated resonances, respec- 
tively 
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FIG. 7: Energy spectrum versus the magnetic field for two 
atoms in separated harmonic traps in the vicinity of a 
Feshbach resonance at Bo = 22.42mT. The resuhs are calcu- 
lated for 7? = 2 , uj = 2-RX lOOkHz, d = 0.0, 0.1, 2.1, 4.0/, with 
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FIG. 8: Energy spectrum versus the magnetic field for two 
^"K atoms in separated harmonic traps in the vicinity of a 
Feshbach resonance at Bo = 22.42mT. The results are calcu- 
lated for T] = 0.7, LP = 27r X lOOkHz, d = 0.0, 0.1, 2.1, 4.0?, 
with I — ^JhJiJiWz)- 



IV. WAVE FUNCTIONS 

A. Analytic expressions 

We turn to the analysis of eigenfunctions. For an 
eigenenergy e the corresponding wave function ^ ^ is given 
by Eq. pO)) . or its integral representation ([72]) . In fact 
those formulas arc not very convenient for the numerical 
evaluation of ^'£(r), because the former involves double 
summation which converges slowly, while the latter is 
valid only at e < 0. For the purpose of numerical calcu- 
lations we derive yet another representation containing 
only a single summation. To this end we express the de- 
nominator in Eq. pi?]) in terms of the integral pT|) . and 
we perform only a single summation either over n or fc 
quantum numbers. Summing over n with the help of the 
generating function (IA3p leads to 



7r3/2 



-r(- 



2-q 



(22) 



where U{a,b,x) denotes confluent hypergcomctric func- 
tion. Both r(a;) and U{a,b,x) are well defined arbitrary 
values of the parameters a, b and x, hence this represen- 
tation is applicable for all values of the energy. One can 
easily check that for small r the wave function is singu- 
lar and behaves as 1/r. which is an expected behavior 
resulting from the application of delta pseudopotcntial. 

When we instead sum over the quantum number k in 
Eq. Uni), with the help of (jAi]) we obtain 



7r3/2 



ti£^coth(0 + ^ 



Vl-e-2* 



(23) 



In this case the integral cannot be expressed in terms of 
a closed formula, except the special case of d = (23| . 



B. Results 



different trap separations. We observe a similar struc- 
ture of avoided crossings as in Figs. S] and [5l For d = 
where only states with s-wave symmetry contribute, the 
spectrum is regular without avoided crossings. Increasing 
d results in the appearance of the oscillator states with 
odd Uz, creating avoided crossing with the rest of the 
states. Finally for sufficiently large d, the trap-induced 
resonances start to appear. 



In FigurcOwe present two example wave functions cal- 
culated for a = /, 7] = 1; in the vicinity of an avoided 
crossing at d = 2.1/. The corresponding eigenvalues, 
E = — 0.048?ia; and E sa 0.327?kj are shown in the panel 
b). Both wave functions are multiplied by r in order 
to obtain a regular behavior at r —^ 0. The wave func- 
tions presented in panel a) and c) are, respectively, an- 
tisymmetric and symmetric combinations of the bound 
and of the trap states. When, moving away from the 
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trap-induced resonance region, the wave functions be- 
come dominated by the single component locahzed ei- 
ther at r = dz or at r = 0, taking the character of trap 
vibrational or the bound state, respectively. 



V. EXAMPLE APPLICATIONS OF THE 
TRAP-INDUCED RESONANCES 

To illustrate applicability of our results we consider 
two example systems, where the trap-induced shape res- 
onances can be applied and can bring some advantages 
over the standard schemes based on collisions in the same 
potential well. First we consider the quantum phase gate, 
involving a pair of separated particles, and controlled by 
external magnetic field. As a second example we present 
a scheme allowing for a quantum state control of the rel- 
ative motion in the trap. In the case of atoms in an 
optical lattice this method allows for a coherent transfer 
of particle into higher Bloch bands. In our analysis we do 
not calculate the full quantum dynamics for the consid- 
ered schemes, but rather focus on the regime where the 
dynamics can be described in the adiabatic or diabatic 
approximation, using the Landau-Zener theory. 



A. Quantum gate with two separate atoms 

We consider a relatively simple quantum gate build up 
of two atoms trapped in the separate trapping potentials. 
The particles are initially prepared in the ground state 
of the external motion in the trap, and the qubit states 
are encoded in the internal spin (hyperfine) states of the 
atoms. Our scheme assumes that the particles remain at 
fixed positions during the gate operation, and the state 
dependent dynamics required for the quantum gate op- 
eration is gained by applying an external magnetic field. 
This differs from the majority of other theoretical pro- 
posals when the particles are moved to perform the gate 
operation, which can reduce the fidelity due to the exci- 
tations to higher vibrational states. 

We denote the qubit states as |0)l {\0)r) and 
(|l)i?) for particles in the left (right) trap, respectively. 
The two-qubit states are denoted as = \i)L\j)R, for 
i,j = 0, 1. The goal is to realize the conditional phase 
gate that is described by the following truth table: 



100) 


^ e'-^'^'lOO) ^ 


|00), 


101) 


— > e'-^nilOl) — > 


101), 


|10) 


— > e'-^ii'llO) — > 


|io), 


111) 


— > e'-^iMll) — > 


e"^|ll). 



(24) 



The first transformation Ub is the actual gate process 
performed with application of an external magnetic field, 
where (jjij describe the phases gained by the two-qubit 
states \i)L\j)R- The second transformation Us represent 



the single qubit operations that can be applied to undo 
the dynamic phases in all, except the single channel [3^ , 
which gains the phase (j) = (f>QQ + (jin — (/)io — 4>oi ^ . For 
(j) = TT a so-called phase gate is achieved, which combined 
with single qubit rotations, form a controUed-NOT. This 
is particularly interesting since controUed-NOT together 
with single qubit operations constitute a universal set for 
quantum computation (40| . 

The main idea behind our quantum gate is to exploit 
the dependence of the avoided crossings positions on the 
individual two-qubit states This dependence stems 
from the fact that each of the two-qubit states, which 
can be treated as separate scattering channel, will exhibit 
distinct structure of the trap-induced shape resonances, 
resulting from differences in the actual positions of the 
molecular states (except states |01) and |10) in the case of 
like atoms). In this way, Feshbach resonances in each of 
the channel state appear at different values of the mag- 
netic field, which allows to realize the state-dependent 
dynamics. In some range of the magnetic field, certain 
two-qubit states can exhibit an avoided crossings (see 
Figlinia), whereas the other states will not experience a 
resonance (see FiglTUlb'). 

By changing the magnetic field we pass the avoided 
crossing adiabatically, starting from large toward the 
small magnetic fields, in order to create molecular states 
in the two-qubit states for which the Feshbach resonance 
occur (see Fig[T0la). Then we reverse the process and 
return with the magnetic field to its initial value. At the 
same time the remaining two-qubit states do not experi- 
ence the resonance (see Fig[TUlb) and they only acquire 
the phase due to the Zeeman shifts of the hyperfine levels. 

For simplicity in our analysis we consider the evolu- 
tion in the adiabatic regime, where the dynamics in the 
vicinity of an avoided crossing can be described in the 
framework of the Landau Zener theory. In the case 
of an adiabatic process 0^ = J^' E,j{B(t))dt, where 
the whole operation takes place between ti and tf and 
Eij{B) ~ Ei{B) + Ej{B) denotes the eigenenergy for a 
two-qubit state 

The main limitation to the gate fidelity results from 
the nonadiabatic transitions to other states of external 
motion in the trap, while passing an avoided crossing. 
This yields a lower bound for the gate operation time. It 
can be estimated by applying the Landau-Zener formula 
for the probability of an adiabatic transfer across the 
avoided crossing 

Pad = l-cxp(^-(^^ (25) 



^ In the case of unlike atoms, the single qubit operations can be 
performed by applying appropriate microwave pulses, that can 
address selectively the atoms in different traps, because of differ- 
ent hyperfine structure of particles. In the case of like species, 
one can use the fact that the atoms are spatially separated, which 
allows to address them selectively [Tsll . 
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FIG. 9: (Color online) Panels a) and c): wave function of two interacting atoms confined in isotropic harmonic potentials 
separated by the distance d — 2.11, for the scattering length a — I, and corresponding eigenergies E = 0.048/iaj and E = 0.327hui, 
respectively. Panel b): energy spectrum versus trap separation d. The red (upper) and the green (lower) dots indicate the 
position of eigenstates shown in panels a) and c), respectively. All the lengths are expressed in harmonic oscillator units 



and assuming the linear ramp of the magnetic field: v = 
\dB/dt\ = const. Here AE denotes the level separation 
at the avoided crossing, dE/dB — dEi/dB — dE2/dB 
where Ei{B), E2{B) are the energies of the two levels 
far from the avoided crossing, and v = \dB/dt\ is a ramp 
speed. For our example calculation of the Feshbach res- 
onance for '"^K atoms (see Fig. [7]), we estimate that for 
trap separations d/l = 2.1 the transfer fidelity F — 0.999 
can be achieved with the speed limited to v = 8mT/s, 
resulting in the total operation around 1ms. The value 
of = 0.999 is the order of minimal fidelity required to 
apply the quantum error correction schemes. The gate 
operation time may be further improved by the optimal 
control techniques [l^, [4l|, H^l . Finally we note that the 
gate operation time will crucially depend on the geome- 
try of the traps, which trivially stems from the fact that 
for large rj states are elongated along the symmetry axis, 
which leads to stronger overlap of the ground-state wave 
functions in the traps and as a consequence to stronger 
avoided crossings. 



B. Transport to higher bands using trap-induced 
shape resonances 

As a second example we consider a scheme for a quan- 
tum state control of a pair of atoms in a harmonic trap. It 
allows for a transfer of particles between arbitrary states 
of the relative motion in a trap. A similar method, based 
on the use of Feshbach resonances has been applied ex- 
perimentally, to transfer a mixture of fermions trapped 
in a deep optical lattice into higher Bloch bands [lj| . In 
our scheme we combine the technique of a magnetic field 
sweep with the use of trap-induced shape resonances, to 




b) 



magnetic field 



magnetic field 



FIG. 10: (Color online) Quantum gate for atoms in separate 
traps controlled by external magnetic fields. Depending on 
the two-qubit states the atoms would experience (panel a) or 
not (panel b) the Feshbach resonance, which leads to the state 
dependent evolution. 



obtain a method for an efficient transfer of atoms into ar- 
bitrary excited state of the relative motion. This can be 
applied, for instance, to probe states of atoms in hi ghe r 
Bloch bands, using the similar setup as in Ref. [ij]. 
This requires that external potential for different spin 
states can be controlled separately, e.g. by means of spin- 
dependent lattice potentials. 

We assume that initially the atoms are prepared in the 
ground state of the trap, and the magnetic field is set at 
the attractive side of a Feshbach resonance. Our scheme 
consists of the following steps (see Fig. [TT|) : a) magnetic 
field sweep across Feshbach resonance resulting in an adi- 
abatic conversion of a pair of atoms into a molecule; b) 
increasing the energy of a molecule by changing the trap 
separation from to some finite value d, before reaching 
the avoided crossing c) conversion of molecule into an ex- 
cited state of atoms in the trap, using the sweep across 
Feshbach resonance d) bringing back the traps to an ini- 



10 



tial positions {d = 0). Steps a) and c) assume constant 
trap separation, whereas b) and d) assume constant value 
of the magnetic field. In principle one could exchange the 
order of steps b) and c), converting a molecule directly 
into excited state of the trap, and later change the value 
of the magnetic field from repulsive to attractive interac- 
tions. Fig. [TT] shows, however, that such a scenario could 
be more difficult to realize, since the strength of avoided 
crossing is of the same order for low-lying and highly ex- 
cited states, when changing the trap distance d instead 
of a magnetic field (compare panels b and c) 

Using the Landau-Zener theory, we made an esti- 
mation for an adiabatic transfer from the vibrational 
ground-state to the third excited state of a relative mo- 
tion. We took parameters for two ''^"K atoms, near a 
Feshbach resonance at Bq = 22.421 niT, assuming spher- 
ically symmetric trap (?; = 1). The main limitation for 
the rate of transfer results from the dynamics during the 
stage c), where the first two avoided crossings have to 
be traversed diabatically, whereas the third one adiabat- 
ically. We investigated the simplest case of a constant 
sweep rate across the Feshbach resonance. For a sweep 
at rate v = \dB/dt\ « 0.2mT/s, the multiple-crossing 
Landau-Zener theory [i^, predicts that atoms at the 
stage c) can be transferred from the molecular to the 
third excited state with efficiency greater than 75%. In 
contrast at stage a) , the avoided crossing is much stronger 
and, for instance, conversion with probability 99% is ob- 
tained with much faster rate v = 0.35T/s. Steps b) and 
d) can be also performed relatively fast, with speed lim- 
ited by the time scale gi ven by t he level separation in the 
trap d <C /w, with I = yJh/JjuJ). Finally all the steps can 
be further optimized applying the optimal control tech- 
niques, which should significantly reduce the time of the 
whole process. For instance a transfer of a particle in a 
harmonic trap to a distance of the order of few harmonic 
oscillator lengths can be performed on the scale of a sin- 
gle trap period T = 27r/a;, by appropriate optimizing the 
transfer process [isj . 



VI. SUMMARY 

In this paper we have presented analytical solutions for 
the system of two atoms trapped in harmonic potentials 
separated by an arbitrary distance d. The atom-atom 
interaction have been modeled with the Fermi pseudopo- 
tential, including the energy-dependent scattering length. 
This accounts for the presence of tight traps and scatter- 
ing resonances. We have investigated the properties of 
the energy spectrum and of the wave functions for differ- 
ent trap separations, scattering lengths, and aspect ratios 
of the cylindrically symmetric trap. The spectrum ex- 
hibit avoided crossings resulting from the resonances be- 
tween the molecular and trap states. For isotropic traps 
(77 = 1) the trap vibrational states at large d, correlate 
at small d with the states of two interacting atoms with 
s and p-wave symmetry. One can explain this behavior 




magnetic field trap separation 




magnetic field trap separation 



FIG. 11: (Color online) Coherent transfer of a pair of atoms 
into an excited state; a) conversion of pair of atoms into 
molecules by a magnetic field sweep; b) transfer of a molec- 
ular state to the region of an avoided crossing (trap-induced 
shape resonance) by changing trap separation; c) conversion 
of a molecule into an excited state in the trap by a magnetic 
field sweep; d) transferring the state back to the regime of a 
zero trap separation. We end with two particles in an excited 
state of a or p-wave symmetry. 



in the framework of the perturbation theory. Finally we 
illustrated our analytical results with two physical ex- 
amples, where the trap-induced shape resonances can be 
applied. The first one was a quantum phase gate for two 
separated particles, based on a control of the magnetic 
field. The second one was a scheme for a coherent trans- 
port of particles into higher Bloch bands. 

In our work we have not considered the possibility 
of the two different trapping frequencies for interacting 
atoms. This can be particularly important for atoms of 
different species or in the case of the ion-atom collisions. 
The inclusion of the two different trapping frequencies 
can lead to the additional resonances between center-of- 
mass and relative motions, that should be present already 
at zero trap separations. 
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APPENDIX A: GENERATING FUNCTIONS FOR 
HARMONIC OSCILLATOR STATES 

The harmonic oscihator functions in one and two- 
dimensional traps are given by. respectively, 



4>„,o(/0,</?) = "-L^XTyp"), m = 0, (Al) 



7r4 V 



(A2) 



In derivation of the integral representation of the wave 
function, we have applied the generating functions for 
Laguerre 

oo 

5^L„(x)z" = (l-z)-iexp(^), (A3) 

n=0 

and Hermite polynomials 



A:=0 



(A4) 



APPENDIX B: DERIVATION OF THE 
RECURRENCE RELATIONS 

Let's consider the difference between values of T(x) 
separated by ?/ 



T(x) - T{x + r])= dt 



_ ,2 l-e-*/2 



(Bl) 



After the substitution y ~ tanh(|;) this can be trans- 
formed into the following form 



T{x)-H^+r^) = 2r// dye-''-yy-''\l-yf'^-\l+y)-^-. 
Ja 

(B2) 

The integral on the r.h.s. can be evaluated analytically 
with the help of [36l | 



= B(i/,A)$i(i/,p,A + i/;/3;-Ai) (B3) 

valid for Re(A) > 0, Re(t/) > 0, |arg(l - /3)| < tt. This 
yields 

Tix)^T{x + r]) = B (i,2a;) $i {\,2x,2x + l\-l,-d'^) , 

(B4) 



Here, $i(a, 6, c, x, y) is the confluent hypergeometric 
function of two variables, which can be defined by the 
series fSol 



^i{a,b,c,x,y) = ^ 



(a),„+„(6) 

X y . 



?n,n— ^ ' 



(B5) 



Evaluation of $i (a, 6, c, x, y) in terms of a double scries is 
not convenient, and in the numerical calculations we have 
applied another result, which can be derived by applying 
the binomial expansion 



{-xfV{k-a) 
k\T{-a) 



(B6) 



to the terms (1 - y)^^"^ or (1 + y)-^^ in (|B2|. In this 
case the remaining integral can be calculated with the 

help of [ail 



dxx''-^{u-xY-^e'^'' 

-- Bill, iy)u^'+''-^iFi{iy; ^i + v- fiu), ^ > 0, > 

(B7) 



In the way described above we obtain two series expan- 
sions for J-{x) — !F{x + Vj): 



T{x) - T[x + On 



(B8) 



with 



r(n + 2^ 
^ \} ' r(2a;)n! 



(2a;)r 

X B(2a::,n-H i)ii^i(n-f i;2x + n+ (B9) 



and 



T{x) - T{x + ri) = 



r(2x; 



n=0 



(BIO) 



with 



For large n, the first series (|B8p converges rather slowly: 
On + dn+i exp(— d^)/?!^ for n odd. In numerical cal- 
culations it is easier to apply the second series (jB10|) . 
which has the asymptotic behavior given by: 6„ 
V2^{n + 2a;)2— 3/2 /2"+4--i/2. 
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